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PREFACE 

This  user's  guide  documents  a  computer  program  that  uses  a  plane  and 
axisymmetric  finite  element  method  of  analysis  for  solution  of  steady-state 
seepage  problems.  This  report  provides  updated  information  on  the  original 
two-dimensional  seepage  program  developed  for  the  U.  S.  Army  Engineer  Division 
Lower  Mississippi  Valley  (LMVD) .  This  report  supersedes  the  original  report 
which  was  published  as  U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES) 
Miscellaneous  Paper  K-73-4,  dated  May  1973,  entitled  "A  Plane  and  Axisymmetric 
Finite  Element  Program  for  Steady-State  and  Transient  Seepage  Problems."  The 
work  in  enhancing  the  program  and  updating  this  report  was  accomplished  as 
part  of  the  applications  support  provided  by  the  WES  Automatic  Data  Processing 
(ADP)  Center  to  LMVD. 

The  work  in  modifying  the  computer  program  and  preparing  the  user's 
guide  was  done  by  Mr.  Fred  T.  Tracy,  Chief,  Research  and  Development  Software 
Group,  ADP  Center. 

Work  on  the  program  was  coordinated  with  LMVD  and  the  U.  S.  Army  Engi¬ 
neer  Districts,  St.  Louis  (LMS)  and  Vicksburg  (LMK) .  Liaison  was  maintained 
between  LMVD,  LMS,  LMK,  and  WES  by  means  of  informal  meetings  and  telephone 
conversations  with  Mr.  Tony  Young,  LMVD;  Mr.  Tom  Wolff,  LMS;  and  Mr.  Eugene 
Wardlaw,  LMS,  all  of  which  provided  technical  guidance  in  defining  the  desired 
program  modifications.  Mr.  Wardlaw  was  also  very  helpful  in  testing  the 
program. 

Dr.  N.  Radhakrishnan,  Special  Technical  Assistant,  ADP  Center,  and 
Mr.  Paul  K.  Senter,  Computer-Aided  Design  Group,  ADP  Center,  coordinated  the 
work. 

Commanders  and  Directors  of  WES  during  this  period  were  COL  Nelson  P. 
Conover,  CE,  and  COL  Tilford  C.  Creel,  CE.  Technical  Director  was  Mr.  F.  R. 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


MultiDl 


To  Obtain 


cubic  feet  per  minute 

0.2831685 

cubic  meters  per  minute 

feet 

0.3048 

meters 

pounds  (mass)  per  cubic  foot 

16.01846 

kilograms  per  cubic  meter 

AXI SYMMETRIC  FINITE  ELEMENT  PROGRAM 


USER'S  GUIDE  FOR  A  PLANE  AND 

FOR  STEADY-STATE  SEEPAGE  PROBLEMS 


PART  I :  INTRODUCTION 


Background 


1.  The  finite  element  method  (FEM)  is  an  outstanding  numerical  tool  for 
solving  two-dimensional  seepage  problems.  It  is  naturally  suited  to  solving 
problems  involving  complex  geometries  or  nonhomo geneous  and  anisotropic  soils. 
Unconfined  flow  problems  can  also  be  solved  by  using  a  FEM  program  to  automat¬ 
ically  determine  the  position  of  the  phreatic  surface. 

2.  To  take  advantage  of  the  FEM  for  Lower  Mississippi  Valley  Division 
(LMVD)  seepage  problems,  a  computer  program  and  report,  "A  Plane  and  Axisymmet 
ric  Finite  Element  Program  for  Steady-State  and  Transient  Seepage  Problems," 
were  written  in  1973.*  The  program  provided  the  following  unique  capabilities 

a.  Allowing  the  free  surface  to  cross  boundaries  of  different  soil 
layers . 

b.  Solving  a  partially  confined  and  partially  unconfined  problem. 

c.  Relieving  restrictions  on  how  to  number  the  grid  for  unconfined 
flow  problems. 

d.  Computing  the  flows  at  the  entrance  and  exit  boundaries. 

3.  However,  some  problems  were  noticed  in  the  program: 

a.  Inputting  incorrect  phreatic  surface  data  caused  the  program  to 
stop  executing. 

b.  Transient  and  steady-state  problems  were  computed  in  the  same 
way.  Thus,  the  steady-state  solution  was  inefficient,  and  most 
of  the  problems  solved  by  LMVD  using  this  program  were  steady- 
state  problems . 

c.  Triangular  elements  were  not  allowed  with  the  free  surface,  and 
the  quadrilateral  elements  had  to  be  numbered  a  certain  way. 

4.  Further,  some  improvements  were  suggested  to  make  the  program  more 
useful : 

a.  Provide  automatic  convergence  for  steady-state  problems. 


*  Miscellaneous  Paper  K-73-4,  U.  S.  Army  Engineer  Waterways  Experiment  Sta 
tion,  Vicksburg,  Miss. 


b.  Simplify  input  boundary  codes  for  unconfined  flow  problems. 

5.  Therefore,  a  modified  version  of  the  original  seepage  program  has 
been  written  for  steady-state  problems  only.  This  new  version  has  the  follow¬ 
ing  advantages: 

a.  Automatic  convergence. 

b.  Easy  determination  of  boundary  codes. 

c.  No  initial  guess  of  the  phreatic  surface  required. 

d.  No  restrictions  on  how  the  grid  is  numbered  or  whether  triangu¬ 
lar  elements  are  used. 


6.  This  report  is  basically  an  update  of  the  report  (Miscellaneous 
Paper  K-73-4)  on  the  original  program.  Thus,  it  is  a  user's  guide  for  the  new 
version  of  the  seepage  program.  Minor  differences  in  input  data  between  the 
old  version  and  this  new,  updated  version  are  given. 
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\  PART  II:  USER'S  GUIDE 


Computer  Program  Requirements 

7.  The  program  is  written  in  FORTRAN  IV.  It  requires  40K  words  of 
core  memory,  four  scratch  disc  files  (logical  units  01,  02,  04,  and  07),  and  a 
permanent  disc  file  for  postprocessor  output  (logical  unit  12) .  Its  Engineer¬ 
ing  Computer  Program  Library  (ECPL)  number  is  704-F3-R0-011 . 

Consistent  Units 

I 

8.  Before  the  actual  data  are  described,  it  should  be  emphasized  that 
data  must  be  input  in  consistent  units.  For  example,  if  head  is  specified  in 
feet  and  time  in  minutes,  then  the  following  units  must  be  used: 

j  a.  Permeability — ft/min. 

b.  Discharge  velocity — ft/min. 

c.  (x,y)  node  point  data--ft. 

Thus,  if  length  (L)  is  initially  given  in  feet  and  time  (T)  in  minutes,  all 
subsequent  input  (and  output)  data  must  use  the  same  units. 

Description  of  Input  Data 

9.  The  following  input  data  cards  are  required  for  each  problem.  Any 
number  of  problems  can  be  solved  by  the  same  computer  run  by  simply  stacking 
problem  data  sets  sequentially.  In  the  section  to  follow,  two  items  for  each 
type  of  data  are  supplied  to  aid  users  who  are  familiar  with  FORTRAN  language 
and  FORMAT  statements.  The  first  item  gives  the  format  specifications  (shown 
in  brackets  after  the  side  headings)  under  which  the  data  are  input.  The  sec¬ 
ond  item  states  the  names  (variables)  used  in  the  program  to  retain  the  data. 
Title  card  [13A6,  A2] 

10.  This  is  an  identification  card  for  the  problem: 

Column  Variable  _ Description _ 

1-80  HED  Alphanumeric  identification 


Any  legal  FORTRAN  character  is  permissible.  The  title  appears  on  the  printed 
output  of  the  program. 


Control  card  [415.  IX.  A4.  F10.0 


This  card  controls  the  size  and  type  of  problem 


Column 


11-15 


16-20 


22-25 


26-35 


Variable 


NUMNP 


NUMEL 


NUMMAT 


NFLCD 


DATUM 


Description 


Number  of  node  points.  NUMNP 
cannot  be  greater  than  1000 

Number  of  elements.  NUMEL 
cannot  be  greater  than  950 

Number  of  soil  types.  NUMMAT 
cannot  be  greater  than  12 

Number  of  discharge  velocity 
cards  (see  paragraph  19) 

Type  of  problem.  Input  PLNE 
for  plane  flow  or  AXSY  for 
axisymmetric  flow 

Elevation  of  datum  (L) .  The 
datum  is  defined  by 

H  =  —  +  Y  -  D 


H  is  the  total  head,  P 
is  the  pressure,  y  is 
the  density  of  water,  Y 
is  the  y  coordinate,  and 
D  is  the  datum.  The  datum 
is  an  arbitrary  constant 
and  is  typically  set  to 
zero.  However,  in  well 
problems  it  can  be  set  to 
the  height  of  the  headwater 
level,  thus  making  the  to¬ 
tal  head  negative  and  cor¬ 
responding  to  drawdown 

This  list  is  shorter  than  that  for  the  old  version  of  the  program  because  the 
transient  data  of  time  increment  (DELTAT)  and  print  interval  (IPRNT)  have  been 
removed.  Also,  a  declaration  of  whether  the  problem  is  confined  or  unconfined 
(LTYPE)  is  no  longer  needed  since  the  program  determines  this  automatically. 
Infiltration  rate  (RINFL)  has  also  been  removed. 

Soil  cards  [15,  2F10.01 

12.  These  cards  describe  the  characteristics  of  soil  types  in  the  prob¬ 
lem.  One  card  is  required  for  each  soil  type: 


Column 


Variable 


Description 


1-5  I  Soil  type  number.  Soil  cards  must 

be  in  ascending  order  with  re¬ 
spect  to  soil  type  number 

6-15  XK1(I)  First  principal  permeability  (L/T) 

16-25  XK2(I)  Second  principal  permeability  (L/T) 

The  only  difference  in  this  list  and  the  old  one  is  that  the  storage  coeffi¬ 
cient  for  each  material  type  (STOR(I))  is  omitted  because  this  information  was 
needed  only  for  transient  problems. 

Nodal  point  cards  [15,  12,  13,  3F10.0] 

13.  These  cards  give  required  information  for  the  nodes: 


Column  Variable 


Description 


1-5  N 

6-7  IN 


8-10  NBC(N) 


11-20  X(N) 

21-30  Y(N) 


Nodal  point  number 

Boundary  condition  flag  for  gener¬ 
ated  nodes  of  a  line  segment,  the 
first  node  of  which  is  N.  If 
IN  =  0,  all  boundary  condition 
information  is  set  to  zero.  If 
IN  =  1,  the  type  of  boundary  con¬ 
dition  for  generated  nodes  is  set 
to  the  value  at  node  N.  The 
boundary  value  for  generated 
nodes  is  determined  by  linear  in¬ 
terpolation  if  IN  =  1  or  is  set 
to  zero  if  JN  =  0 

Boundary  condition  flag.  NBC  =  0 
for  an  internal  node  or  an  imper¬ 
vious  boundary  node;  NBC  =  1  for 
a  specified  head  on  a  boundary; 
and  NBC  =  -1  for  a  specified  flow 
on  a  boundary.  For  unconfined 
flow  problems,  use  NBC  =  2  for 
the  exit  face  (see  Figure  1). 

Be  sure  to  include  all  nodes 
from  the  top  of  the  face,  down 
to  and  including  the  node  at  the 
tailwater  level 

x  value  (L)  of  node  point 

y  value  (L)  of  node  point 


(Continued) 
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Column 


Variable 


Description 


31-40  FX(N)  Boundary  value.  FX  represents  head 

(L)  for  positive  NBC  and  flow  per 
2 

unit  length  (L  /T)  for  plane 
problems  or  flow  per  unit  radian 
3 

(L  /T)  for  axisymmetric  problems 
for  negative  NBC 


NBC  =  0 


Figure  1.  Unconfined  flow  boundary  codes 


One  card  per  node  is  required  unless  a  set  of  consecutively  numbered  nodes  is 
located  at  equal  intervals  along  a  straight-line  segment  (see  Figure  2).  In 
this  case,  nodal  point  data  for  the  interior  nodes 
of  the  line  segment  (nodes  5-8)  are  generated  if 
node  cards  for  these  points  are  omitted.  Thus,  for 
the  example  in  Figure  2,  the  card  for  node  4  should 
be  followed  by  a  card  for  node  9.  The  coordinates 
of  the  generated  nodes  are  found  by  linearly  inter-  Figure  2.  Line  of  nodes 
polating  between  the  coordinates  of  the  first  node 

(4)  and  the  last  node  (9)  of  the  line  segment.  The  type  of  boundary  condition 
for  the  generated  nodes  is  set  either  to  zero  or  to  the  value  of  the  first 
node  (see  variable  IN  above).  Boundary  values  for  generated  nodes  either  are 
determined  by  linearly  interpolating  between  the  first  and  last  nodes  of  a 
line  segment  or  are  set  to  zero. 
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14.  Node  cards  must  be  numbered  in  ascending  order  with  respect  to  node 
number,  and  the  number  of  the  first  node  must  be  1.  A  data  card  for  both  the 
first  and  last  nodes  must  always  be  supplied.  The  old  version  of  the  program 
requested  a  rate  of  change  of  boundary  condition  value  for  each  node 
(DELFX(N)).  This  item  has  been  dropped  from  the  list  because  it  is  used  only 
for  transient  problems.  The  meanings  of  the  NBC  values  have  also  been 
simplified. 

Phreatic  surface  data 

15.  The  old  version  of  the  program  required  an  initial  guess  as  to  the 
phreatic  surface.  The  new  version  determines  this  automatically,  so  these 
data  are  no  longer  required. 

Element  cards  [615,  F10.0] 

16.  These  cards  contain  the  information  about  the  elements: 


Column 

Variable 

Description 

1-5 

J 

Element  number 

6-10 

NP(1,J) 

First  node  of  element  J.  The  nodes 
must  be  in  counterclockwise  order 
as  shown  in  Figure  3.  In  this 
example,  the  nodes  could  be  ren¬ 
dered  5,  10,  11,  6 

11-15 

NP(2,J) 

Second  node  of  element  J 

16-20 

NP(3, J) 

Third  node  of  element  J 

21-25 

NP(4,J) 

Fourth  node  of  element  J 

26-30 

NP(5,J) 

Soil  type  number 

31-40 

ANG(J) 

Angle  (in  degrees)  between  the 

first  principal  permeability  and 
the  x  axis  as  shown  in  Figure  4; 
k  and  k2  represent  the  two  prin' 


cipal  permeabilities 


Figure  3.  Counterclock¬ 
wise  numbering 


Figure  4.  Principal 
permeabilities 
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One  card  per  element  is  needed  unless  a  group  of  elements  satisfies  the  follow¬ 
ing  conditions: 

a.  The  soil  type  and  two  directions  of  the  principal  permeabili¬ 
ties  are  the  same  for  all  elements  of  a  group. 

b.  Element  numbers  of  a  group  are  consecutive. 

c.  Node  numbers  for  elements  N+l  of  a  group  can  be  generated  by 
incrementing  the  node  numbers  of  element  N  by  1 . 


In  this  case,  element  data  for  all  but  the  first  element  of  a  group  are  gener¬ 
ated  if  cards  for  these  elements  are  omitted.  For  example,  elements  3-5  of 
the  configuration  in  Figure  5  form  a  group.  The  data  cards  for  element  3 
should  therefore  be  followed  by  a  data  card 

for  element  6.  The  soil  type  and  directions  4  - - 1'2 

of  the  principal  permeabilities  for  generated 

5  0 

elements  are  set  to  the  respective  values  of 

the  first  element  (element  3)  of  a  group.  3  - - -  n 

Node  numbers  for  generated  elements  are  pro- 

4  7 

duced  by  successively  incrementing  the  respec¬ 
tive  node  numbers  of  the  first  element  of  a  2  - - -  io 

group  by  1. 

3  6 

17.  Element  cards  must  be  numbered  in 

ascending  order  with  respect  to  element  num-  t  - - 1 9 

5 

ber,  and  the  number  of  the  first  element  must 

be  1.  A  data  card  for  both  the  first  and  Figure  5.  Sample  grid 


last  elements  must  always  be  supplied. 

18.  During  initial  setup  of  the  problem,  the  user  should  attempt  to 
assign  node  numbers  to  the  elements  in  such  a  way  as  to  minimize  the  differ¬ 
ence  between  the  largest  and  the  smallest  node  numbers  of  each  element.  This 
so-called  "greatest  difference"  cannot  be  greater  than  59.  The  old  version  of 
the  program  required  special  numbering  of  elements ,  and  no  triangular  elements 
were  allowed  wherever  the  phreatic  surface  would  potentially  go.  In  the  new 
version,  these  restrictions  no  longer  apply. 

Discharge  velocity  cards  1215,  F10.01 


19.  The  discharge  velocity  cards  are  used  to  specify  discharge  velocity 
along  straight-line  segments  connecting  two  nodes: 
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Column 


Variable 


Description 


1-5  K  First  node  of  line  segment 

6-10  L  Second  node  of  line  segment 

11-20  FLRT  Normal  component  of  discharge  ve¬ 

locity  (L/T)  on  line  segment.  A 
positive  value  for  FLRT  indicates 
that  the  flow  is  entering  the 
system;  a  negative  value  indi¬ 
cates  exit 

One  card  for  each  line  segment  is  required.  The  total  number  of  discharge  ve¬ 
locity  cards  must  equal  NFLCD.  If  NFLCD  is  zero,  no  cards  are  needed.  NBC(K) 
and  NBC(L)  must  be  set  to  -1  in  the  nodal  point  data. 

Description  of  Output  Data 

20.  Output  is  in  the  form  of  printed  information  and  a  postprocessor 
file.  Four  groups  of  printed  information  are  supplied: 

a.  Input  data. 

b.  Heads  and  flows  at  nodes. 

c.  Discharge  velocities  (flow  rates)  for  the  elements. 

d.  Phreatic  surface  and  iteration  information  for  unconfined  flow 
problems . 

Input  data 

21.  The  problem  identification,  soil  information,  nodal  point  informa¬ 
tion,  element  information,  and  flow  rate  data  are  first  printed.  These  data 
are  included  for  all  generated  nodes  and  elements.  A  printed  copy  of  the  in¬ 
put  data  is  valuable  for  verifying  that  the  intended  data  were  supplied. 

Heads  and  flows 

22.  The  head  at  every  node  is  printed.  Since  seepage  results  from  a 
head  differential  rather  than  from  the  magnitude  of  head,  the  percentage  of 
net  head  is  also  printed  for  every  node.  Flow  entering  or  leaving  the  flow 
region  is  provided  for  every  boundary  node.  As  shown  in  Figure  6,  flow  at  a 
node  is  the  quantity  of  water  per  unit  time  per  unit  length  (plane  flow)  or 
per  unit  radian  (axisymmetric  flow)  crossing  the  boundary  through  the  effec¬ 
tive  area  surrounding  the  node.  A  positive  flow  indicates  flow  entering  the 
system,  and  a  negative  flow  indicates  exit. 


Figure  6.  Flow 


Discharge  velocities 

23.  The  discharge  velocities  at  Xe,  Ye  in  the  and  k^  directions,  re- 

spectively,  are  printed  for  each  element  in  the  flow  region.  The  value  X  ,  Y 

is  the  average  of  the  x's  and  y's,  respectively,  of  the  element  corner  nodes. 

If  the  phreatic  surface  crosses  a  group  of  elements  as  shown  in  Figure  7,  the 

coordinates  of  the  nodes  of  these  elements  are  temporarily  changed  to  coincide 

with  the  phreatic  surface.  Elements  outside  the  flow  region  are  therefore 
-0  -0 

eliminated,  and  X  ,  Y  is  altered  for  some  elements. 


MODIFIED  GRID 


Figure  7.  Phreatic  surface  intersecting  the  grid 


24.  The  magnitude  and  direction  of  the  discharge  velocity  are  also 
printed.  The  direction  is  specified  by  the  counterclockwise  angle  between  the 
x  axis  and  the  discharge  velocity  direction. 

Phreatic  surface  and 
iteration  information 

25.  Phreatic  surface.  The  status  of  each  node  with  respect  to  the 


PART  III:  SAMPLE  PROBLEMS 


Problem  1:  Unconfined  Flow  in  a  Bank  with  Vertical  Sides 


After  Sudden  Drawdown — Dupuit's  Problem 

Problem  description 

28.  The  problem  deals  with  a  bank  having  vertical  sides  where  the  head¬ 
water  and  tailwater  are  initially  the  same  until  the  tailwater  level  is  sud¬ 
denly  lowered  (see  Figure  8).  The  problem  consists  of  determining  the  phre¬ 
atic  surface,  head  distribution,  and  quantity  of  flow  after  steady-state 
conditions  have  been  reached. 


Figure  8.  Dupuit's  problem 


Values  of  the  constants  are 


H  =  100  ft* 

h  =  20  ft 
o 

L  =  100  ft 
k  =  0.1  ft/min 


Finite  element  grids 


29.  Two  different  FE  grids  are  used  to  solve  this  problem:  a  small 
rectangular  grid  (Figure  9),  and  a  large  triangular  grid  (Figure  10). 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  3. 
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Figure  10.  Large  triangular  grid 
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Input  cards 

30.  A  listing  of  the  input  cards  for  the  rectangular  grid  is  given 


below. 


DUPUIT'S  PROBLEM 


36 

25 

1 

0 

PLNE 

0. 

1 

.1 

.1 

1  1 

1 

0. 

0. 

100 

6 

1 

0. 

100. 

100 

7 

20. 

0. 

12 

20. 

100. 

13 

40. 

0. 

18 

40. 

100. 

19 

60. 

0. 

24 

60. 

100. 

25 

80. 

0. 

30 

80. 

100. 

31 

1 

100. 

0. 

20 

32  1 

2 

100. 

20. 

20 

36 

2 

100. 

100. 

100 

1 

1 

7 

8 

2 

1 

6 

7 

13 

14 

8 

1 

11 

13 

19 

**0 

14 

1 

16 

19 

25 

26 

20 

1 

21 

25 

31 

32 

26 

1 

25 

29 

35 

36 

30 

1 

Output 

31.  A  listing  of  the  output  for  the  rectangular  grid  is  presented  below. 

PLANE  FLOW  PROBLEM 


DUPUIT'S  PROBLEM 


NUMBER  OF  NODAL  POINTS - 36 

NUMBER  OF  ELEMENTS - 25 

NUMBER  OF  DIFF.  MATERIALS—  1 

ELEVATION  OF  DATUM -  0, 

MATERIAL  PROPERTIES 
MAT  K1  K2 

1  0.100E  00  0 . 100E  00 
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NODE  POINT  INFORMATION 


NODE 

BC 

X 

Y 

FLOW-HEAD 

1 

1 

0. 

0. 

100.00 

2 

1 

0* 

20.00 

100.00 

3 

1 

0. 

40.00 

100.00 

4 

1 

0. 

60.00 

100.00 

5 

1 

0. 

80.00 

100.00 

6 

1 

0. 

100.00 

100.00 

7 

0 

20 » 00 

0. 

0. 

8 

0 

20 » 00 

20.00 

0. 

9 

0 

20.00 

40.00 

0. 

10 

0 

20.00 

60.00 

0. 

11 

0 

20.00 

80.00 

0. 

12 

0 

20.00 

100.00 

0. 

13 

0 

40.00 

0. 

0. 

14 

0 

40.00 

20.00 

0. 

15 

0 

40.00 

40.00 

0. 

16 

0 

40.00 

60.00 

0. 

17 

0 

40.00 

80.00 

0. 

18 

0 

40.00 

100.00 

0. 

19 

0 

60.00 

0. 

0. 

20 

0 

60.00 

20.00 

0. 

21 

0 

60.00 

40.00 

0* 

22 

0 

60.00 

60.00 

0. 

23 

0 

60.00 

80.00 

0. 

24 

0 

60.00 

100.00 

0. 

25 

0 

80.00 

0. 

0. 

26 

0 

80.00 

20.00 

0. 

27 

0 

80.00 

40.00 

0. 

28 

0 

80.00 

60.00 

0. 

29 

0 

80.00 

80.00 

0. 

30 

0 

80.00 

100.00 

0. 

31 

1 

100.00 

0. 

20.00 

32 

2 

100.00 

20.00 

20.00 

33 

2 

100.00 

40.00 

40.00 

34 

2 

100.00 

60.00 

60.00 

35 

2 

100.00 

80.00 

80.00 

36 

2 

100.00 

100.00 

100.00 

NODAL  FLOWS  AND  HEADS  POSITION  OF  PHREATIC  SURFACE 
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0.2000E  02  0.  %  -0.U92E  01  * 

0.2000E  02  0.  X  -0.2938E  01  * 

0.4000E  02  25.0  X  -0.6832E  00  * 

0.4035E  02  25.4  X  *  100.00  40.35 


ELEMENT  FLOW  RATES 


ELMT 

VI 

V2 

P-AXIS  ANG 

RES  V 

DIR  OF  V 

1 

0.605E-01 

-0.864E-03 

0. 

0.605E-01 

-0 .818E 

00 

2 

0 « 572E-01 

-0 . 242E-02 

0. 

0 . 573E-01 

-0.242E 

01 

3 

0 • 512E-01 

-0  *  362E-02 

0. 

0.513E-01 

-0.404E 

01 

4 

0 • 429E-01 

-0 ♦ 463E-02 

0. 

0 » 432E-01 

-0.615E 

01 

5 

0 . 336E-01 

-0.469E-02 

0. 

0  *  339E-01 

-0.795E 

01 

6 

0.640E-01 

-0.280E-02 

0. 

0.641E-01 

-0.250E 

01 

7 

0  » 600E-01 

-0. 778E-02 

0. 

0  >  605E-01 

-0.738E 

01 

8 

0.531E-01 

-0.112E-01 

0. 

0.543E-01 

-0.119E 

02 

9 

0.448E-01 

-0 « 136E-01 

0. 

0.469E-01 

-0.168E 

02 

10 

0 • 401E-01 

-0 • 121E-01 

0. 

0 ♦ 419E-01 

-0.168E 

02 

11 

0.719E-01 

-0.543E-02 

0. 

0 . 721E-01 

-0.432E 

01 

12 

0.661E-01 

-0« 150E-01 

0. 

0 . 678E-01 

-0.128E 

02 

13 

0.561E-01 

-0.209E-01 

0. 

0.599E-01 

-0.204E 

02 

14 

0 « 468E-01 

-0.214E-01 

0. 

0 . 515E-01 

-0.246E 

02 

15 

0 • 454E-01 

-0 « 140E-01 

0. 

0  *  475E-01 

-0.1 72E 

02 

16 

0.871E-01 

-0 • 120E-01 

0. 

0. 879E-01 

-0.782E 

01 

17 

0 « 752E-01 

-0.262E-01 

0. 

0*796E-01 

-0.192E 

02 

18 

0.614E-01 

-0.335E-01 

0. 

0.700E-01 

-0.286E 

02 

19 

0 • 547E-01 

-0.263E-01 

0. 

0.607E-01 

-0.257E 

02 

21 

0.116E  00 

-0 • 84AE-02 

0. 

0.117E  00 

-0.416E 

01 

22 

0.914E-01 

-0.665E-01 

0. 

0.113E  00 

-0.360E 

02 

23 

0 • 569E-01 

-0 . 416E-01 

0. 

0.705E-01 

-0.362E 

02 

ITERATION 

NO. 

= 

1 

DELTA 

H 

MAX. 

= 

0* 10000000E 

31 

ITERATION 

NO. 

S 

2 

DELTA 

H 

MAX. 

X 

0.57387257E 

00 

ITERATION 

NO. 

s 

3 

DELTA 

H 

MAX. 

= 

0.20409732E 

02 

ITERATION 

NO. 

s 

4 

DELTA 

H 

MAX. 

X 

0 » 12305222E 

01 

ITERATION 

NO. 

s 

5 

DELTA 

H 

MAX. 

X 

0 . 1 1060864E 

02 

ITERATION 

NO. 

X 

6 

DELTA 

H 

MAX. 

= 

0.63306427E 

00 

ITERATION 

NO. 

ae 

7 

DELTA 

H 

MAX. 

X 

0.93912921E 

01 

ITERATION 

NO. 

s 

8 

DELTA 

H 

MAX. 

X 

0.47176266E 

00 

ITERATION 

NO. 

ss 

9 

DELTA 

H 

MAX. 

X 

0 . 271 17305E 

01 

ITERATION 

NO. 

s 

10 

DELTA 

H 

MAX. 

X 

0.77009296E 

00 

ITERATION 

NO. 

X 

11 

DELTA 

H 

MAX* 

= 

0.29948235E 

00 

ITERATION 

NO. 

X 

12 

DELTA 

H 

MAX. 

X 

0. 10039425E 

00 

ITERATION 

NO. 

X 

13 

DELTA 

H 

MAX. 

X 

0. 18034458E 

00 

ITERATION 

NO. 

X 

14 

DELTA 

H 

MAX. 

X 

0.55501938E- 

-01 

29  4  24 


M  M  144  124 


Figure  11.  Modified  grid  plot 


A  -  9.3444E  42 
I  •  9.4444E  92 

C  •  0.5444E  42 
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F  •  9.8944E  92 
G  •  9.8444E  92 
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Figure  12.  Equipotential  line  plot 
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Figure  15.  Equipotential  line  plot 
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Figure  16.  Velocity  vector  plot 


k,  =  1  X  10-3  CM/SEC 
k2  =  1/9  x  10-3  CM/SEC 


Figure  17.  Weir  problem 


Figure  19.  Model  of  sheet  pile 

Input  cards 

36.  A  listing  of  the  input  data  is  given  below.  Note  that  the  grid 
being  generated  has  a  card  for  each  node  and  a  card  for  each  element.  Note 
also  that  units  of  feet  and  minutes  are  used  in  the  problem. 


A 

t 


r 

A7 

0 

-94.000 

-25.600 

0. 

0, 

i 

48 

0 

-70.000 

-6.500 

0. 

0, 

£ 

49 

0 

-48.000 

9.600 

0. 

0. 

50 

0 

-28.000 

23.000 

0. 

0. 

51 

0 

-10.000 

34.000 

0. 

0. 

a 

. 

52 

0 

6.000 

42.900 

0. 

0. 

53 

0 

20.000 

50.000 

0. 

0, 

h 

54 

0 

32.000 

55.600 

0. 

0. 

\ 

55 

1 

42.000 

60.000 

70.000 

0. 

56 

0 

-120.000 

-60.000 

70.000 

0. 

57 

0 

-94.000 

-36.300 

0. 

0. 

i 

58 

0 

-70.000 

-16.000 

0. 

0. 

59 

0 

-48.000 

1.200 

0. 

0. 

> 

60 

0 

-28.000 

15.600 

0. 

0. 

» 

61 

0 

-10.000 

27.500 

0. 

0. 

1 

62 

0 

6.000 

37.200 

0. 

0. 

1 

63 

0 

20.000 

45.000 

0. 

0. 

64 

0 

32.000 

51.200 

0. 

0. 

65 

0 

42.000 

56.100 

0. 

0. 

66 

1 

50.000 

60.000 

70.000 

0. 

• 

67 

0 

-94.000 

-47.000 

0. 

0, 

. 

68 

0 

-70.000 

-25.500 

0. 

0. 

69 

0 

-48.000 

-7.200 

0. 

0. 

70 

0 

-28.000 

8.200 

0. 

0. 

k 

71 

0 

-10.000 

21.000 

0. 

0. 

w 

72 

0 

6.000 

31.500 

0. 

0. 

u 

c 

73 

0 

20.000 

40.000 

0. 

0. 

k 

74 

0 

32.000 

46.800 

0. 

0. 

A 

75 

0 

42.000 

52.200 

0. 

0. 

76 

0 

50.000 

56.500 

0. 

0. 

• 

77 

0 

55.000 

60.000 

0, 

0. 

78 

0 

-70.000 

-35.000 

0, 

0. 

-’ 

79 

0 

-48.000 

-15.600 

0. 

0. 

BO 

0 

-28.000 

0.800 

0. 

0. 

B1 

0 

-10.000 

14.500 

0. 

0. 

»■ 

82 

0 

6.000 

25.800 

0. 

0. 

•\ 

63 

0 

20.000 

35.000 

0. 

0. 

*  - 
k  • 

84 

0 

32.000 

42.400 

0. 

0. 

85 

0 

42.000 

48.300 

0. 

0. 

► 

86 

0 

50.000 

53.000 

0. 

0. 

L 

87 

0 

55.000 

56.750 

0. 

0. 

88 

0 

59.999 

60.000 

0. 

0. 

b'. 

89 

0 

-48.000 

-24.000 

0. 

0. 

»  %‘ 
h*’ 

90 

0 

-28.000 

-6.600 

0. 

0. 

91 

0 

-10.000 

8.000 

0. 

0. 

s 

t 

r 

92 

0 

6.000 

20.100 

0. 

0. 

93 

0 

20.000 

30.000 

0. 

0. 

94 

0 

32.000 

38.000 

0. 

0. 

95 

0 

42.000 

44.400 

0. 

0. 

96 

0 

50.000 

49.500 

0. 

0. 

97 

0 

55.000 

53.500 

0. 

0. 

* 


29 


49 

0 

65.000 

43.500 

0. 

0 

150 

0 

70.000 

47.500 

0. 

0 

.51 

0 

75.000 

50.000 

0. 

0 

[52 

0 

7C.000 

52.500 

0. 

0 

L53 

0 

50.000 

25.000 

0. 

0 

L54 

0 

55.000 

30.750 

0. 

0 

155 

0 

60.000 

36.000 

0. 

0 

156 

0 

65.000 

40.750 

0. 

0 

L57 

0 

70.000 

45.000 

0. 

0 

L  56 

0 

75.000 

47.500 

0. 

0 

L59 

0 

75.000 

52.500 

0. 

0 

160 

0 

80.000 

50.000 

0. 

0 

L61 

0 

70.001 

55.000 

0. 

0 

L62 

0 

55.000 

27.500 

0. 

0 

163 

0 

60.000 

33.000 

0. 

0 

L64 

0 

65.000 

38.000 

0. 

0 

L65 

0 

70.000 

42.500 

0. 

0 

L  66 

0 

75.000 

45.000 

0. 

0 

167 

0 

80.000 

47.500 

0. 

0 

6  S 

0 

75.000 

55.000 

0. 

0 

L  69 

0 

80.000 

52.500 

0. 

0 

0 

85.000 

50.000 

0. 

0 

L71 

0 

70.001 

57.500 

0. 

0 

172 

0 

60.000 

30.000 

0. 

0 

L73 

0 

65.000 

35.250 

0. 

0 

L74 

0 

70.000 

40.000 

0. 

0 

L75 

0 

75.000 

42.500 

0. 

0 

[76 

0 

80.000 

45.000 

0. 

0 

[77 

0 

85.000 

47.500 

0. 

0 

[78 

0 

75.001 

57.500 

0. 

0 

[79 

0 

80.000 

55.000 

0. 

0 

0 

85.000 

52.500 

0. 

0 

[81 

0 

90.000 

50.000 

0. 

0 

.82 

1 

70.001 

60.000 

60.000 

0 

.83 

0 

65.000 

32.500 

0. 

0 

[84 

0 

70.000 

37.500 

0. 

0 

185 

0 

75.000 

40.000 

0. 

0 

[86 

0 

80.000 

42.500 

0. 

0 

87 

0 

85.000 

45.000 

0. 

0 

38 

0 

90.000 

47,500 

0. 

0 

[89 

1 

75.001 

60.000 

60.000 

0 

90 

0 

80.001 

57.500 

0. 

0 

91 

0 

85.000 

55.000 

0. 

0 

[92 

0 

90.000 

52.500 

0. 

0 

L93 

0 

95.000 

50.000 

0. 

0 

94 

0 

70.000 

35.000 

0. 

0 

[95 

0 

75.000 

37.500 

0. 

0 

[96 

0 

80.000 

40.000 

0. 

0 

L97 

0 

85.000 

42.500 

0. 

0 

98 

0 

90.000 

45.000 

0. 

0 

199 

0 

95.000 

47.500 

0. 

0 

200 

1 

80 

.001 

60 

.000 

60.000 

201 

0 

85 

.001 

57 

.500 

0. 

202 

0 

90 

.000 

55 

.000 

0. 

203 

0 

95 

.000 

52 

.500 

0. 

204 

0 

100 

.000 

50 

.000 

0. 

205 

1 

85 

.001 

60 

.000 

60.000 

206 

0 

90 

,000 

57 

.500 

0. 

207 

0 

95 

.000 

55 

,000 

0. 

208 

0 

100 

.000 

52 

.500 

0. 

209 

0 

105 

.000 

52 

.500 

0. 

210 

1 

90 

.001 

60 

.000 

60.000 

>11 

0 

95 

.000 

57 

.500 

0. 

212 

0 

100 

.000 

55 

.000 

0. 

213 

0 

105 

.000 

55 

.000 

0. 

214 

0 

110 

.000 

55 

.000 

0. 

>15 

1 

95 

.000 

60 

.000 

60.000 

>16 

0 

100 

.000 

57 

.500 

0. 

>17 

0 

105 

.000 

57 

.500 

0. 

218 

0 

110 

.000 

57 

,500 

0. 

*19 

0 

115 

.000 

57 

.500 

0. 

*20 

1 

100 

.000 

60 

.000 

60.000 

*21 

1 

105 

.000 

60 

.000 

60.000 

*22 

1 

110 

.000 

60 

.000 

60.000 

>23 

1 

115 

.000 

60 

.000 

60.000 

!24 

1 

120 

.000 

60 

.000 

60.000 

1 

66 

76 

87 

77 

1 

26.565 

T 

A 

76 

86 

97 

87 

1 

26.565 

3 

86 

96 

107 

97 

1 

26.565 

4 

96 

106 

117 

107 

1 

26.565 

5 

106 

116 

126 

117 

1 

26.565 

6 

116 

125 

134 

126 

1 

26,565 

7 

125 

133 

141 

134 

1 

26.565 

8 

133 

140 

147 

141 

1 

26.565 

o 

140 

146 

154 

147 

1 

26.565 

10 

146 

153 

162 

154 

1 

26.565 

11 

77 

87 

98 

88 

1 

26,565 

12 

87 

97 

108 

98 

1 

26.565 

13 

97 

107 

118 

108 

1 

26.565 

14 

107 

117 

127 

118 

1 

26.565 

15 

117 

126 

135 

127 

1 

26.565 

16 

126 

134 

142 

135 

1 

26.565 

17 

134 

141 

148 

142 

1 

26.565 

18 

141 

147 

155 

148 

1 

26.565 

19 

147 

154 

163 

155 

1 

26.565 

20 

154 

162 

172 

163 

1 

26.565 

21 

88 

98 

109 

99 

1 

26.565 

22 

98 

108 

119 

109 

1 

26.565 

23 

108 

118 

128 

119 

1 

26.565 

24 

118 

127 

136 

128 

1 

26.565 

25 

127 

135 

143 

136 

1 

26.565 

26 

135 

142 

149 

143 

1 

26.565 

L-, 

'  ‘ 

.  *.  -  .  •  .  -  -  -V  * '  T. i"  i-:  i-'  *"  i i  i  j  v  .  w  .• 

y 

ft 

,  ^ 

27 

142 

148 

156 

149 

i 

26.565 

28 

149 

155 

164 

156 

i 

26.565 

29 

155 

163 

173 

164 

i 

26.565 

30 

163 

172 

183 

173 

i 

26.565 

31 

99 

109 

120 

110 

i 

26.565 

32 

109 

119 

129 

120 

i 

26.565 

■ 

33 

119 

128 

137 

129 

i 

26.565 

34 

128 

136 

144 

137 

i 

26.565 

“*  *■ 

35 

136 

143 

150 

144 

i 

26.565 

g 

36 

143 

149 

157 

150 

i 

26.565 

1 

37 

149 

156 

165 

157 

i 

26.565 

38 

156 

164 

174 

165 

i 

26.565 

39 

164 

173 

184 

174 

i 

26.565 

40 

173 

183 

194 

184 

i 

26.565 

/• 

41 

153 

146 

139 

145 

i 

26.565 

jr 

42 

146 

140 

132 

139 

i 

26.565 

1 

43 

140 

133 

124 

132 

i 

26.565 

44 

133 

125 

115 

124 

i 

26.565 

w* 

45 

125 

116 

105 

115 

i 

26.565 

46 

116 

106 

95 

105 

i 

26.565 

47 

106 

96 

85 

95 

i 

26.565 

48 

96 

86 

75 

85 

i 

26.565 

BT 

49 

86 

76 

65 

75 

i 

26.565 

• 

50 

76 

66 

55 

65 

i 

26.565 

V, 

51 

145 

139 

131 

138 

i 

26.565 

52 

139 

132 

123 

131 

i 

26.565 

53 

132 

124 

114 

123 

i 

26.565 

g 

54 

124 

115 

104 

114 

i 

26.565 

* 

55 

115 

105 

94 

104 

i 

26.565 

56 

105 

95 

84 

94 

i 

26.565 

»  * 

57 

95 

85 

74 

84 

i 

26.565 

i;v 

58 

85 

75 

64 

74 

i 

26.565 

.  ’ 

59 

75 

65 

54 

64 

i 

26.565 

L 

60 

65 

55 

45 

54 

i 

26.565 

61 

138 

131 

122 

130 

i 

26.565 

*  *.• 

62 

131 

123 

113 

122 

i 

26.565 

63 

123 

114 

103 

113 

i 

26.565 

r_ 

64 

114 

104 

93 

103 

i 

26.565 

65 

104 

94 

83 

93 

i 

26.565 

.C 

66 

94 

84 

73 

83 

i 

26.565 

It 

67 

84 

74 

63 

73 

i 

26.565 

68 

74 

64 

53 

63 

i 

26.565 

69 

64 

54 

44 

53 

i 

26.565 

70 

54 

45 

36 

44 

i 

26.565 

71 

130 

122 

112 

121 

i 

26.565 

•» 

72 

122 

113 

102 

112 

i 

26.565 

* 

73 

113 

103 

92 

102 

i 

26.565 

74 

103 

93 

82 

92 

i 

26.565 

75 

93 

83 

72 

82 

i 

26.565 

76 

83 

73 

62 

72 

i 

26.565 

77 

73 

63 

52 

62 

i 

26.565 

. 
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78  63 

79  53 

80  44 

81  121 
82  112 

83  102 

84  92 

85  82 

86  72 

87  62 

88  52 

89  43 

90  35 

91  111 

92  101 

93  91 

94  81 

95  71 

96  61 

97  51 

98  42 

99  34 

100  27 

101  100 

102  90 

103  80 

104  70 

105  60 

106  50 

107  41 

108  33 

109  26 

110  20 

111  89 

112  79 

113  69 

114  59 

115  49 

116  40 

117  32 

118  25 

119  19 

120  14 

121  78 

122  68 

123  58 

124  48 

125  39 

126  31 

127  24 

128  18 


26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 


129 

13 

9 

5 

8 

130 

9 

6 

3 

5 

131 

67 

57 

46 

56 

132 

57 

47 

37 

46 

133 

47 

38 

29 

37 

134 

38 

30 

22 

29 

135 

30 

23 

16 

22 

136 

23 

17 

11 

16 

137 

17 

12 

7 

11 

138 

12 

8 

4 

7 

139 

8 

5 

2 

4 

140 

5 

3 

1 

2 

141 

220 

216 

217 

221 

142 

216 

212 

213 

217 

143 

212 

208 

209 

213 

144 

208 

204 

209 

209 

145 

221 

217 

218 

222 

146 

217 

213 

214 

218 

147 

213 

209 

214 

214 

148 

222 

218 

219 

223 

149 

218 

214 

219 

219 

150 

223 

219 

224 

224 

151 

182 

171 

178 

189 

152 

171 

161 

168 

178 

153 

161 

152 

159 

168 

154 

152 

144 

151 

159 

155 

189 

178 

190 

200 

156 

178 

168 

179 

190 

157 

168 

159 

169 

179 

158 

159 

151 

160 

169 

159 

200 

190 

201 

205 

160 

190 

179 

191 

201 

161 

179 

169 

180 

191 

162 

169 

160 

170 

180 

163 

205 

201 

206 

210 

164 

201 

191 

202 

206 

165 

191 

180 

192 

202 

166 

180 

170 

181 

192 

167 

210 

206 

211 

215 

168 

206 

202 

207 

211 

169 

202 

192 

203 

207 

170 

192 

181 

193 

203 

171 

215 

211 

216 

220 

172 

211 

207 

212 

216 

173 

207 

203 

208 

212 

174 

203 

193 

204 

208 

175 

144 

150 

158 

151 

176 

150 

157 

166 

158 

177 

157 

165 

175 

166 

178 

165 

174 

185 

175 

179 

174 

184 

195 

185 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 


26 

.565 

26 

.565 

26 

.565 

26 

.565 

26 

.565 

26 

.565 

26 

.565 

26 

.565 

26 

.565 

26 

.565 

26 

.565 

26 

.565 

26 

.565 

26 

.565 

26 

.565 

26 

.565 

26 

.565 

26 

.565 

26 

.565 

26 

,565 

26, 

,565 

26, 

565 

26, 

,565 

26, 

,565 

26. 

565 

26. 

565 

26. 

565 

26. 

565 

26. 

565 

26. 

565 

26. 

565 

26. 

565 

26. 

565 

26. 

565 

26. 

565 

26, 

565 

26, 

565 

26, 

565 

26. 

565 
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=\ n  «v  -  --tv 

180 

184 

194 

195 

195 

181 

151 

158 

167 

160 

182 

158 

166 

176 

167 

183 

166 

175 

186 

176 

184 

175 

185 

196 

186 

185 

185 

195 

196 

196 

186 

160 

167 

177 

170 

187 

167 

176 

187 

177 

188 

176 

186 

197 

187 

189 

186 

196 

197 

197 

190 

170 

177 

188 

181 

191 

177 

187 

198 

188 

192 

187 

197 

198 

198 

193 

181 

188 

199 

193 

194 

188 

198 

199 

199 

195 

193 

199 

204 

204 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 

26.565 


Output 


37.  A  listing  of  the  output  is  now  given. 


PLANE  FLOW  PROBLEM 


WEIR  PROBLEM 


NUMBER  OF  NODAL  POINTS - 224 

NUMBER  OF  ELEMENTS - - 195 

NUMBER  OF  DIFF  .  MATERIALS—  1 
ELEVATION  OF  DATUM— .  0. 


MATERIAL  PROPERTIES 
MAT  K1  K2 

1  0.197E-02  0  » ? 19E-03 


37 


Jl 


NODE  POINT  INFORMATION 


NODE  BC 


FLOW-HEAD 


1 

1 

-120.00 

60.00 

70.00 

0 

-120.00 

48.00 

70.00 

3 

1 

-94.00 

60.00 

70.00 

4 

0 

-120.00 

36.00 

70.00 

5 

0 

-94.00 

49.30 

0. 

6 

1 

-70.00 

60.00 

70.00 

7 

0 

-120.00 

24.00 

70.00 

8 

0 

-94.00 

38.60 

0. 

9 

0 

-70.00 

50.50 

0. 

10 

1 

-48.00 

60.00 

70.00 

11 

0 

-120.00 

12.00 

70.00 

12 

0 

-94.00 

27.90 

0. 

13 

0 

-70.00 

41.00 

0. 

14 

0 

-48.00 

51.60 

0. 

15 

1 

-28.00 

60.00 

70.00 

16 

0 

-120.00 

0. 

70.00 

17 

0 

-94.00 

17.20 

0. 

18 

0 

-70.00 

31.50 

0. 

19 

0 

-48.00 

43.20 

0. 

20 

0 

-28.00 

52.60 

0. 

21 

1 

-10.00 

60.00 

70.00 

22 

0 

-120.00 

-12.00 

70.00 

23 

0 

-94.00 

6.50 

0. 

24 

0 

-70.00 

22.00 

0. 

25 

0 

-48.00 

34.80 

0. 

26 

0 

-28.00 

45.20 

0. 

27 

0 

-10.00 

53.50 

0. 

28 

1 

6.00 

60.00 

70.00 

29 

0 

-120.00 

-24.00 

70.00 

30 

0 

-94.00 

-4.20 

0. 

31 

0 

-70.00 

12.50 

0. 

32 

0 

-48.00 

26.40 

0. 

33 

0 

-28.00 

37.80 

0. 

34 

0 

-10.00 

47.00 

0. 

35 

0 

6.00 

54,30 

0. 

36 

1 

20.00 

60.00 

70.00 

37 

0 

-120.00 

-36.00 

70.00 

146 

$ 

50.00 

28.50 

0. 

147 

0 

55  00 

34.00 

0. 

143 

0 

£0.00 

39,00 

0. 

149 

0 

£5.00 

43.50 

0. 

150 

0 

70 . 00 

47.50 

0. 

151 

0 

75.00 

50,00 

A. 

,j 

152 

0 

70.00 

52,50 

0 . 

153 

0 

50.00 

25.00 

0 . 

* 

154 

a 

55.00 

30.75 

0. 

155 

0 

60.00 

36.00 

0. 

156 

0 

65.00 

40.75 

0. 

■H 

157 

A 

70.00 

45.00 

0. 

158 

0 

75.00 

47,50 

A  . 

15? 

.■v 

1  1 

V 

75.00 

52.50 

0. 

*« 

160 

0 

80.00 

50.00 

0. 

161 

0 

70.00 

55,00 

0. 

162 

0 

55.00 

27.50 

0. 

•* 

163 

0 

60.00 

33.00 

0. 

i 

164 

0 

65.00 

38 . 00 

0. 

\ 

165 

0 

70.00 

42,50 

A 

V  « 

\ 

166 

0 

75.00 

45 . 00 

0. 

167 

A 

90.00 

47.50 

0. 

168 

o 

75.00 

55  00 

0, 

fa 

16? 

A 

30.00 

52.50 

0. 

1 

170 

0 

35.00 

50.00 

0. 

* 

171 

A 

70.00 

57.50 

0. 

172 

0 

60.00 

30,00 

0. 

7 

* 

173 

0 

65.00 

35.25 

0. 

■ 

174 

A 

70.00 

40.00 

0. 

i 

175 

0 

75.00 

42.50 

0. 

* 

176 

0 

80.00 

45.00 

0. 

177 

0 

85.00 

47.50 

0. 

J 

178 

0 

75.00 

57.50 

0. 

17? 

0 

30.00 

55.00 

0. 

180 

0 

85.00 

52.50 

0, 

i 

181 

0 

90.00 

50.00 

0. 

182 

1 

70.00 

60.00 

60,00 

■ 

183 

0 

65.00 

22.50 

0, 

* 

194 

0 

70.00 

37.50 

0, 

185 

0 

75.00 

40,00 

0. 

186 

0 

90.00 

42.50 

0. 

187 

A 

85.00 

45.00 

0, 

188 

A 

90.00 

47.50 

0, 

18? 

1 

75.00 

60,00 

60.00 

190 

A 

80.00 

57.50 

0. 

■ 

1?1 

0 

85.00 

55.00 

0. 

■i 

« 

192 

0 

90.00 

52.50 

0, 

S 

193 

0 

95.00 

50,00 

0. 

j 

194 

0 

70.00 

35.00 

0. 

. 

195 

A 

75.00 

37,50 

0. 

196 

0 

90,00 

40,00 

0. 

i 


\ 


► 


r.*.- 


I 


197 

0 

85.00 

42.50 

0. 

198 

0 

90.00 

45.00 

0. 

199 

0 

95.00 

47.50 

0. 

200 

1 

80.00 

60.00 

60.00 

201 

0 

85.00 

57.50 

0. 

202 

0 

90.00 

55.00 

0. 

203 

0 

95.00 

52.50 

0. 

204 

0 

100.00 

50.00 

0. 

205 

1 

85.00 

60.00 

60.00 

206 

0 

90.00 

57.50 

0. 

207 

0 

95.00 

55.00 

0. 

208 

0 

100.00 

52.50 

0. 

209 

0 

105.00 

52.50 

0. 

210 

1 

90.00 

60.00 

60.00 

211 

0 

95.00 

57.50 

0. 

212 

0 

100.00 

55.00 

0. 

213 

0 

105.00 

55.00 

0. 

214 

0 

110.00 

55.00 

0. 

215 

1 

95.00 

60.00 

60.00 

216 

0 

100.00 

57.50 

0. 

217 

0 

105.00 

57.50 

0. 

218 

0 

110.00 

57.50 

0. 

219 

0 

115.00 

57.50 

0. 

220 

1 

100.00 

60.00 

60.00 

221 

1 

105.00 

60.00 

60.00 

222 

1 

110.00 

60.00 

60.00 

223 

1 

115.00 

60.00 

60.00 

224 

1 

120.00 

60.00 

60.00 

ELEMENT  INFORMATION 


ELMT 

11 

•2 

*3 

1 

66 

76 

87 

2 

76 

86 

97 

3 

86 

96 

107 

4 

96 

106 

117 

5 

106 

116 

126 

6 

116 

125 

134 

7 

125 

133 

141 

8 

133 

140 

147 

9 

140 

146 

154 

#4 

MAT 

ANGLE 

77 

1 

26.6 

87 

1 

26.6 

97 

1 

26.6 

107 

1 

26.6 

117 

1 

26.6 

126 

1 

26.6 

134 

1 

26.6 

141 

1 

26.6 

14-» 

1 

26.6 

42 


w  r  j 


63 

205 

201 

206 

210 

64 

201 

191 

202 

206 

65 

191 

180 

192 

202 

66 

180 

170 

181 

192 

67 

210 

206 

211 

215 

68 

206 

202 

207 

211 

69 

202 

192 

203 

207 

70 

192 

181 

193 

203 

71 

215 

211 

216 

270 

72 

211 

207 

212 

216 

73 

207 

203 

208 

212 

74 

203 

193 

204 

208 

75 

144 

150 

158 

151 

76 

150 

157 

166 

158 

PERCENTAGE  OF 


NODE 

HEAD 

AVAILABLE  HEAD 

FL  Ol' 

1 

0.7000E  02 

100.0  X 

0.9356E-05 

2 

0.6997E  02 

99.7  7. 

3 

0.7000E  02 

100.0  X 

0.5691E-04 

4 

0.6989E  02 

98.9  7. 

5 

0.6993E  02 

99.3  7. 

6 

0.7000E  02 

100.0  X 

0 . 1094E-03 

■» 

t 

0.6974E  02 

97.4  7. 

8 

0.6981E  02 

98.1  7. 

9 

0.6990E  02 

99.0  X 

10 

0.7000E  02 

100.0  7. 

0.1542E-03 

11 

0.6954E  02 

95.4  7. 

12 

0.6964E  02 

96.4  X 

13 

0.6975E  02 

97.5  X 

14 

0.6987E  02 

98.7  % 

15 

0.7000E  02 

100.0  X 

0 . 1955E-03 

16 

0.6928E  02 

92.8  X 

17 

0.6943E  02 

94.3  X 

18 

0.6956E  02 

95.6  X 

19 

0.6970E  02 

97.0  X 

20 

0.6985E  02 

98,5  X 

21 

0.7000E  02 

100.0  X 

0.239.3E-03 

22 

0.6900E  02 

90.0  X 

23 

0.6918E  02 

91.8  X 

24 

0.6933E  02 

93.3  X 

25 

0.6948E  02 

94.8  % 

26 

0.6964E  02 

96.4  X 

27 

0.6981E  02 

98.1  X 

28 

0.7000E  02 

100.0  X 

0.2921E-03 

29 

0.6874E  02 

87.4  X 

30 

0.6891E  02 

89.1  X 

31 

0.6907E  02 

90.7  X 

32 

0.6922E  02 

92.2  X 

33 

0.6939E  02 

93.9  X 

34 

0.6957E  02 

95,7  X 

35 

0.6977E  02 

97.7  % 

36 

0.7000E  02 

100.0  X 

0.3613E-03 

37 

0.6853E  02 

85.3  7. 

38 

0.6867E  02 

86.7  X 

39 

0.6880E  02 

88,0  7. 

40 

0.6894E  02 

89.4  % 

41 

0.6909E  02 

90.9  X 

42 

0.6927E  02 

92.7  7. 

43 

0.6947E  02 

94.7  X 

44 

0.6971E  02 

97.1  X 

45 

0.7000E  02 

100.0  X 

0.4591E-03 

46 

0.6840E  02 

84.0  X 

47 

0.6848E  02 

84,8  X 

48 

0.6856E  02 

85,6  X 

i 


47 


0.6364E 

02 

86.4 

% 

0.6376E 

02 

87.6 

X 

0.6890E 

02 

89,0 

y 

0.6909E 

0  ! 

90.9 

V. 

0.6931E 

02 

93.1 

7. 

0.6961E 

02 

96.1 

7. 

0.7000E 

02 

100.0 

X 

0.6836E 

02 

83.6 

X 

0.6837E 

02 

83.7 

X 

0.6837E 

02 

83.7 

X 

0.6838E 

02 

83.8 

y 

tm 

0.6843E 

02 

84.3 

X 

0.6851E 

02 

35.1 

X 

0.6863E 

02 

86 . 3 

X 

0.6881E 

02 

88.1 

X 

0.6905E 

02 

90.5 

X 

0.6940E 

02 

94.0 

X 

0.7000E 

02 

100.0 

X 

0.6834E 

02 

83.4 

y 

•m 

0.6826E 

02 

82.6 

X 

0.6819E 

02 

81.9 

X 

0.6814E 

02 

81.4 

X 

0.6813E 

02 

81.3 

X 

0.6816E 

02 

81.6 

X 

0.6823E 

02 

82.3 

X 

0.6833E 

02 

83.3 

X 

0.6B43E 

02 

84,3 

X 

0.6845E 

02 

84.5 

•/ 

0.6850E 

02 

85.0 

% 

0.6824E 

02 

82.4 

X 

0.6808E 

02 

80.8 

X 

0.6793E 

02 

79.3 

y 

0.6782E 

02 

78.2 

z 

0.6774E 

02 

77.4 

V 

•9 

0.6770E 

02 

77.0 

X 

0.6768E 

02 

76.8 

X 

0.6767E 

02 

76.7 

X 

0.6766E 

02 

76.6 

X 

0.6777E 

02 

77.7 

y 

•• 

0.6768E 

02 

76.8 

X 

0.6807E 

02 

80.7 

X 

0.6782E 

02 

78.2 

X 

0.6759E 

02 

75.9 

X 

0.6740E 

02 

74.0 

X 

0.6724E 

02 

72 . 4 

V 

9 

0.6712E 

02 

71.2 

X 

0.6703E 

02 

70,3 

X 

0.6700E 

02 

70.0 

X 

0.6714E 

02 

71.4 

X 

0  >  6716E 

02 

71.6 

X 

0.6704E 

02 

70.4 

X 

0 . 6351E-03 


0. 1127E-02 


i 


? 

% 

-* 

* 

* 

r. 

i*. 

i 

i 

'"rf 

k 

l 

>.* 


100 

0.6782E  02 

78.2 

101 

0.6748E  02 

74.3 

102 

0.6717E  02 

71.7 

103 

0.6689E  02 

68.9 

104 

0.6666E  02 

66.6 

105 

0.6648E  02 

64.8 

106 

0.6638E  02 

63.8 

107 

0.6653E  02 

65.3 

108 

0.6664E  02 

66.4 

109 

0.6667E  02 

66.7 

110 

0.6663E  02 

66.3 

111 

0.6749E  02 

74.9 

112 

0.6706E  02 

70.6 

113 

0.6667E  02 

66.7 

114 

0.6631E  02 

63.1 

115 

0.6602E  02 

60.2 

116 

0.6582E  02 

58.2 

117 

0.6591E  02 

59.1 

118 

0.6603E  02 

60.3 

119 

0.6617E  02 

61.7 

120 

0.6628E  02 

62.8 

121 

0.6708E  02 

70.8 

122 

0.6657E  02 

65.7 

123 

0.6610E  02 

61.0 

124 

0.6570E  02 

57.0 

125 

0.6537E  02 

53.7 

126 

0.6532E  02 

53.2 

127 

0.6534E  02 

53.4 

128 

0.6547E  02 

54.7 

129 

0.6571E  02 

57.1 

130 

0.6659E  02 

65.9 

131 

0.6602E  02 

60.2 

132 

0.6552E  02 

55.2 

133 

0.6509E  02 

50,9 

134 

0.6489E  02 

48.9 

135 

0.6470E  02 

47.0 

136 

0.6461E  02 

46.1 

137 

0.6491E  02 

49.1 

138 

0.6606E  02 

60.6 

139 

0.6546E  02 

54.6 

0.6496E  02 

49,6 

141 

0.6465E  02 

46.5 

142 

0.6432E  CC 

43.2 

143 

0.6386E  02 

38.6 

144 

0.6303E  02 

30.3 

145 

0.6550E  02 

55,0 

146 

0.6493E  02 

49.3 

147 

0.6456E  02 

45,6 

148 

0.6418E  02 

41,8 

149 

0.6373E  02 

37.3 

150 

0.6309E  02 

30,9 

% 

% 

X 

7. 

7. 

X 

7. 

7. 

7. 

X 

v 

/• 

X 

7. 

X 

X 

X 

7 


X 

X 

X 

7 

X 

X 

y. 

X 

X 

X 

X 

X 

•/ 

tm 

X 

X 

X 

X 

7. : 

x 

x 

x 

X 

X 

X 

X 

X 

1 

•• 

X 


49 


■V^T 


;  {  .'J 


0.6228E 

02 

22.8 

X 

0.6127E 

02 

12.7 

X 

0.6499E 

02 

49.9 

X 

0.6456E 

02 

45. 6 

X 

0.6414E 

02 

41.4 

X 

0.6368E 

02 

36.8 

X 

0.6313E 

02 

31.3 

X 

0.6250E 

02 

25.0 

X 

0.6190E 

02 

19.0 

X 

0.6190E 

02 

19.0 

X 

0.6052E 

02 

5.2 

X 

0.6463E 

02 

46.3 

X 

0.6416E 

02 

41.6 

X 

0.6370E 

02 

37.0 

X 

0.6318E 

02 

31.8 

X 

0.6265E 

02 

26.5 

X 

0.6211E 

02 

21.1 

X 

0.6096E 

02 

9.6 

z 

0.6159E 

02 

15.9 

X 

0.6160E 

02 

16.0 

X 

0.6012E 

02 

1.2 

X 

0.6424E 

02 

42.4 

X 

0.6374E 

02 

37.4 

X 

0.6324E 

02 

32.4 

X 

0.6275E 

02 

27.5 

X 

0.6227E 

02 

22.7 

X 

0.6180E 

02 

18.0 

X 

0.6034E 

02 

3.4 

X 

0.6115E 

02 

11.5 

X 

0.6135E 

02 

13.5 

X 

0.6135E 

02 

13.5 

X 

0.6000E 

02 

0. 

X 

0.6382E 

02 

38.2 

X 

0.6331E 

02 

33.1 

X 

0.628SE 

02 

28.5 

X 

0.6239E 

02 

23.9 

X 

0.6194E 

02 

19.4 

X 

0.6151E 

02 

15.1 

X 

0.6000E 

02 

0. 

X 

0.6051E 

02 

5.1 

X 

0.6100E 

02 

10.0 

X 

0.6113E 

02 

11.3 

X 

0  •  61 1  IE 

02 

11.1 

X 

0.6339E 

02 

33.9 

X 

0.6294E 

02 

29,4 

X 

0.6248E 

02 

24.8 

X 

0.6204E 

02 

20.4 

X 

0.6161E 

02 

16.1 

X 

0.6121E 

02 

12.1 

X 

0.6000E 

02 

0. 

X 

0.6054E 

02 

5.4 

X 

-0.3593E-04 


-0.2296E-03 


-0.4520E-03 


202 

0.6085E  02 

3,5  2 

203 

0.6094E  02 

9.4  Z 

204 

0.6085E  02 

8.5  % 

205 

0.6000E  02 

0.  7. 

-0.5883E-03 

206 

0.6047E  02 

4-7  7. 

207 

0.6071E  02 

7.1  % 

208 

0.6075E  02 

7.5  7. 

209 

0.6054E  02 

5.4  % 

210 

0.6000E  02 

0.  7. 

-0.5865E-03 

211 

0.6040E  02 

4.0  7. 

212 

0.6058E  02 

5.8  Z 

213 

0.6045E  02 

4.5  7 

214 

0.6028E  02 

2.8  % 

215 

0.6000E  02 

0,  Z 

-0.5019E-03 

216 

0.6033E  02 

3.3  % 

217 

0.6026E  02 

2.6  Z 

218 

0.6019E  02 

i.9 

219 

0.6009E  02 

0.9  Z 

220 

0.6000E  02 

0.  % 

-0.4216E-03 

221 

0 .6000E  02 

0.  % 

-0.3453E-03 

222 

0.6000E  02 

0.  Z 

-0.7668E-0 3 

223 

0.6000E  02 

0.  Z 

-0 , 1787E-03 

224 

0.6000E  02 

0.  z 

-0 . 3207E-04 

ELEMENT  FlOU  RATES 


ELMT 

VI 

V2 

F-AXIS  ANG 

RES  V 

DIR  OF  V 

1 

0. 102E-03 

-0.8B1E-04 

0.266E 

02 

0 1 135E-03 

-0.143E 

02 

2 

0.520E-04 

-0.546E-04 

0.266E 

02 

0.754E-04 

-0.199E 

02 

3 

0.493E-04 

-0.4B5E-04 

0.266E 

02 

0.692E-04 

-0.180E 

02 

4 

0.608E-04 

-0.482E-04 

0.266E 

02 

0.776E-04 

-0.118E 

02 

5 

0.896E-04 

-0.467E-04 

0.266E 

02 

0 . 101E-03 

-0.974E 

00 

6 

0.122E-03 

-0.387E-04 

0.266E 

02 

0. 128E-03 

0.90  IE 

01 

7 

0.139E-03 

-0.264E-04 

0.266E 

02 

0.142E-03 

0.158E 

02 

8 

0 . 138E-03 

-0.159E-04 

0.266E 

02 

0 i 139E-03 

0.200E 

02 

9 

0. 133E-03 

-0.856E-05 

0.266E 

02 

0  * 134E-03 

0.229E 

02 

10 

0.129E-03 

-0.272E-05 

0.266E 

02 

0 , 129E-03 

0.254E 

02 

11 

0.865E-04 

-0.532E-04 

0.266E 

02 

0 . 102E-03 

-0.501E 

01 

12 

0.565E-04 

-0 . 484E-04 

0.266E 

02 

0.744E-04 

-0,140E 

02 

13 

0.477E-04 

-0.500E-04 

0.266E 

02 

0.691E-04 

-0.198E 

02 

14 

0 . 667E-04 

-0.552E-04 

0.266E 

02 

0.866E-04 

-0.130E 

02 

51 


0. 113E-03 

-0.543E-04 

0.2&6E 

02 

0.125E-03 

0.864E 

00 

0.157E-03 

-0.40AE-04 

0.26AE 

02 

0.162E-03 

0.120E 

02 

0.165E-03 

-0.236E-04 

0.266E 

02 

0.16&E-03 

0  * 184E 

0? 

0.153E-03 

-0 . 137E-04 

0.266E 

02 

0 . 153E-03 

0.215E 

02 

0.144E-03 

-0.733E-05 

0.2A&E 

02 

0 . 145E-03 

0.237E 

02 

0.139E-03 

-0.233E-05 

0.266E 

02 

0. 139E-03 

0.256E 

02 

0.690E-04 

-0.417E-04 

0.26AE 

02 

0.806E-04 

-0.460E 

01 

0.376E-04 

-0  *  455E-04 

0.266E 

02 

0.590E-04 

-0.239E 

02 

0.329E-04 

-0.572E-04 

0.266E 

02 

0.660E-04 

-0.335E 

02 

0.738E-04 

-0.701E-04 

0.266E 

02 

0  * 102E-03 

-0.170E 

02 

0.160E-03 

-0.678E-04 

0.266E 

02 

0 . 174E-03 

0.361E 

01 

0.216E-03 

-0.342E-04 

0.266E 

02 

0.219E-03 

0.176E 

02 

0.183E-03 

-0.178E-04 

0.266E 

02 

0.184E-03 

0.210E 

02 

0.165E-03 

-0. 104E-04 

0 . 266E 

02 

0.165E-03 

0.230E 

02 

0.154E-03 

-0.5A1E-05 

0.266E 

02 

0 . 154E-03 

0.245E 

02 

0.148E-03 

-0. 181E-05 

0.266E 

02 

0 . 148E-03 

0.259E 

02 

0.26AE-04 

-0.353E-04 

0.266E 

02 

G.442E-04 

-0.2A4E 

02 

0.276E-05 

-0.497E-04 

0.2&6E 

02 

0.497E-04 

-0.666E 

02 

0.792E-05 

-0.689E-04 

0.266E 

02 

0.693E-04 

-0.700E 

02 

0.771E-04 

-0. 132E-03 

0.266E 

02 

0 » 153E-03 

-0.332E 

02 

0.352E-03 

-0.512E-04 

0.266E 

02 

0.356E-03 

0.1 83E 

02 

0.233E-03 

-0.178E-04 

0.266E 

02 

0.234E-03 

0.222E 

02 

0* 195E-03 

-0. 104E-04 

0.266E 

02 

0.195E-03 

0.235E 

0? 

0. 173E-03 

-0.&26E-05 

0.266E 

02 

0.173E-03 

0.245E 

02 

0.161E-03 

-0 .353E-05 

0.266E 

02 

0.161E-03 

0.253E 

02 

0* 155E-03 

-0. 114E-05 

0.266E 

02 

0. 155E-03 

0.261E 

02 

0. 116E-03 

-0.301E-05 

0.266E 

02 

0.116E-03 

0.251E 

02 

0.118E-03 

-0.940E-05 

0.2A6E 

02 

0.119E-03 

0.220E 

02 

0.119E-03 

-0.170E-04 

0.266E 

02 

0.120E-03 

0.184E 

02 

0.113E-03 

-0.263E-04 

0.266E 

02 

0. 116E-03 

0.134E 

02 

0.947E-04 

-0 . 351E-04 

0.266E 

02 

0.101E-03 

0.619E 

01 

0.701E-04 

-0.409E-04 

0.266E 

02 

0.812E-04 

-0.371E 

01 

0.490E-04 

-0.441E-04 

0.26&E 

02 

0.659E-04 

-0.154E 

02 

0.332E-04 

-0.485E-04 

0.2A6E 

02 

0.588E-04 

-0.290E 

02 

0.111E-04 

-0.590E-04 

0.2&6E 

02 

0.600E-04 

-0.528E 

02 

0.138E-03 

-0.634E-04 

0.266E 

02 

0. 152E-03 

-0.129E 

03 

0.985E-04 

-0.310E-05 

0.266E 

02 

0.985E-04 

0.248E 

02 

0.990E-04 

-0.962E-05 

0.266E 

02 

0.995E-04 

0.210E 

02 

0.957E-04 

-0.168E-04 

0.266E 

02 

0.971E-04 

0.166E 

02 

0.858E-04 

-0.244E-04 

0.26AE 

02 

0.892E-04 

0.107E 

02 

0.689E-04 

-0.309E-04 

0.266E 

02 

0.755E-04 

0.240E 

01 

0.484E-04 

-0.357E-04 

0.266E 

02 

0.601E-04 

-0.983E 

01 

0.262E-04 

-0.395E-04 

0.266E 

02 

0.474E-04 

-0.298E 

02 

0.425E-05 

-0.431E-04 

0.266E 

02 

0.433E-04 

-0.691E 

02 

0.653E-04 

-0.415E-04 

0.266E 

02 

0.774E-04 

-0.121E 

03 

0.818E-04 

-0.247E-04 

0.266E 

02 

0.855E-04 

-0.137E 

03 

0.797E-04 

-0.295E-05 

0.266E 

02 

0.797E-04 

0.244E 

02 

0.789E-04 

-0.908E-05 

0.26&E 

02 

0.794E-04 

0.200E 

02 

0.735E-04 

-0.154E-04 

0.266E 

02 

0.751E-04 

0.147E 

02 

0.626E-04 

-0.215E-04 

0.266E 

02 

0.662E-04 

0.762E 

01 

0.465E-04 

-0.265E-04 

0.266E 

02 

0.535E-04 

-0.310E 

01 

i 


I 

-- 


% 


66 

0.263E-04 

-0 . 301E-04 

0.266E 

02 

0.399E-04 

-0 ,223E 

02 

6? 

0 . 915E-06 

-0.319E-04 

0  -  266E 

02 

0.3ioE-04 

-0.618E 

02 

68 

-0 . 314E-04 

-0.300E-04 

0.266E 

02 

C.434E-04 

-0,1 10E 

03 

69 

-0.490E-04 

-0.222E-04 

0.266E 

02 

0.537E-04 

-0. 129E 

02 

70 

-0.530E-04 

-0.148E-04 

0.266E 

02 

0.550E-04 

-0.138F 

03 

71 

0 . 617E-04 

-0.263E-05 

0.266E 

02 

0 , 618E-04 

0.241E 

02 

72 

0.602E-04 

-0.804E-05 

0.266E 

02 

0 . 608E-04 

0.190E 

02 

73 

0.541E-04 

-0.134E-04 

0.266E 

02 

0.558E-04 

0.127E 
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38.  Figure  20  is  the  equipotential  plot,  and  Figure  21  is  the  velocity 
vector  plot. 

Analysis 

39.  This  problem  can  be  solved  graphically  by  transformation  and  then 
construction  of  a  flow  net.  Using  a  shape  factor  of  the  transformed  problem 

-3 

of  1/2,  the  rate  of  seepage  per  linear  foot  of  structure  is  3.2808  (10  )  cfm. 

-3 

The  flow  from  the  computer  run  is  3.6387  (10  )  cfm,  giving  a  percent  difference 

of  10.3. 
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Figure  20.  Equipotential  line  plot 


APPENDIX  A:  PROCEDURE  FOR  UNCONFINED  FLOW  PROBLEMS 


1.  Determination  of  the  phreatic  surface  in  this  computer  program  is 
based  on  a  modification  of  work  by  Bathe  and  Khashgoftaar.*  An  outline  of  the 
procedure  used  in  the  program  is  as  follows: 

a.  Assume  an  initial  position  of  the  phreatic  surface  at  the  high* 
est  possible  elevation. 

b.  Compute  heads  and  flows  for  each  node,  satisfying  only  the  flow 
=  0  condition  on  the  phreatic  surface. 

c.  Assume  all  points  whose  computed  pressure  is  less  than  zero  to 
be  above  the  phreatic  surface.  To  model  the  "no  flow"  condition 
of  these  nodes,  assume  very  small  permeabilities  in  the  region 
of  these  points. 


d. 


Using  the  new  permeabilities,  compute  a  new  set  of  flows  at  the 

nodes.  All  internal  nodes  should  have  zero  flows  (where  there 

are  no  sources  or  sinks).  Compute  a  set  of  Ah's  to  make  these 

flows  vanish,  and  compute  head  at  node  i  as  h.  =  h.  , . 

+  Ah  .  i  new  i  old 

i 

Repeat  the  iterative  process  of  steps  £  and  d  until  the  maximum 
h  is  of  acceptable  smallness. 


2.  Now,  an  element  whose  nodes  have  all  negative  pressures  is  com¬ 
pletely  above  the  phreatic  surface,  and  an  element  whose  nodes  have  some  posi¬ 
tive  pressures  or  negative  pressures  has  the  phreatic  surface  passing  through 
it.  Linear  interpolation  is  used  to  determine  the  modified  shapes  of  the  af¬ 
fected  elements  (see  Figures  7,  11,  and  14  in  the  main  text)  and  the  position 
of  the  phreatic  surface. 


*  Bathe,  K.-J.,  and  Khashgoftaar,  M.  R.  1979.  "Finite  Element  Free  Surface 
Analysis  without  Mesh  Iteration,"  International  Journal  of  Numerical  and 
Analytical  Methods  in  Geomechanics,  Vol  3,  pp  13-22. 
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